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ABSTRACT 


Transmission  measurements  using  six  wavelengths  were  made 
through  the  edges  of  plumes  for  solid  propellant  rocket  motors 
using  various  propellants  and  motor  geometries.  The  average 
values  obtained  for  the  aluminum  oxide  particles  were  a  Sauter 
mean  diameter  of  0.30+/-. 02  microns,  an  index  of  refraction  of 
1.64+/-. 04  and  standard  deviation  of  1.52+/-. 12  for  the  assumed 
monomodal,  log-normal  size  distribution.  The  results  indicated 
that  the  small  aluminum  oxide  particles  in  the  plume  edge  were 
gamma-Al203,  independent  of  propellant  composition,  motor 
operating  conditions  and  nozzle  geometry.  The  good  correlation 
of  the  data  indicated  that  the  small  particles  can  be  adequately 
represented  by  a  monomodal,  log-normal  distribution. 
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I .  INTRODUCTION 


Aluminum  is  one  of  several  metals  that  can  be  added  to 
solid  rocket  propellants  to  increase  the  specific  impulse. 
The  addition  of  metals  to  the  propellant,  however,  increases 
the  exhaust  signature  of  the  plume.  The  particles  in  the 
plume  are,  in  general,  a  result  of  the  complete  combustion  of 
the  aluminum  [Ref.  1,2].  The  particle  size  distribution 
throughout  the  exhaust  nozzle  and  plume  affects  motor 
performance,  plume  visibility,  infrared  signature,  and  the 
environment  [Ref.  3,4].  The  particle  size  distribution  and 
particle  optical  properties  vary  throughout  the  plume.  Plume 
flowfield  models,  such  as  SPF,  and  plume  radiation  models, 
like  SIRRM,  are  used  to  predict  the  plume  signature,  however 
none  of  the  models  have  been  adequately  validated  [Ref.  5,6]. 

Two  commonly  used  methods  to  analyze  particles  within 
the  plume  are  sample  withdrawal,  and  in  situ  optical 
measurements.  In  the  sample  withdrawal  procedure,  particles 
are  collected  on  a  filter  or  solid  surface,  or  they  are 
removed  through  a  tube  for  analysis.  The  optical  techniques 
measure  the  transmission  or  the  scattering  of  light  that  is 
directed  through  the  plume.  The  later  techniques  have  the 
advantage  of  providing  continuous  measurements  without 
disturbing  the  flowfield. [Ref .  7] 
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Highly  aluminized  propellants  produce  plumes  containing 
large  amounts  of  aluminum  oxide.  Particles  in  the  outer  edges 
of  the  plume  can  then  dominate  the  plume  signature.  Small 
particles  (micron,  submicron)  dominate  in  this  outer  plume 
region  because  the  larger  particles  cannot  turn  with  the  gas 
flow  through  the  nozzle  throat  and  follow  the  wall  profile. 
Thus,  the  optical  properties  and  size  of  the  small  particles 
are  important  for  plume  signature  prediction. [Ref .  2] 

Diagnostic  techniques  used  for  determining  particle  size 
and  optical  properties  generally  have  limited  dynamic  ranges. 
Techniques  that  work  for  submicron  particles  generally  do  not 
work  for  large  particles,  and  visa-versa.  One  technique  that 
has  been  successfully  used  for  submicron  particles  is 
multiple-wavelength  light  transmission  measurements  [Ref. 
3,7].  These  measurements  are  sometimes  augmented  by  light 
scattering  measurements  [Ref.  4,8]. 

In  a  recent  investigation  at  the  Naval  Postgraduate  School, 
the  multiple-wavelength  technique  was  used  to  measure  the 
particles  in  the  edges  of  the  plume  of  a  small  solid  propellant 
rocket  motor.  It  was  found  that  the  particles  in  this  region  were 
dominated  by  small  particles  which  had  an  index  of  refraction  of 
1.63+8%,  and  a  Sauter  mean  diameter  of  0.15m+4%.  Recommended 
improvements  to  that  investigation  were  to  significantly  increase 
the  rate  of  data  acquisition  and  incorporate  additional 
wavelengths,  to  reduce  the  effects  of  time-variation  in  motor 
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operating  conditions,  and  to  provide  more  data  for  improved 
accxiracy  of  tne  technique.  [Ref .  9] 

Since  the  small  particles  in  the  edges  of  the  plume  can 
dominate  the  plume  signature,  it  is  desirable  to  know  whether  or 
not  their  size  and  optical  properties  are  somewhat  "universal”  or 
whether  they  vary  significantly  with  propellant  composition,  motor 
geometry,  and  /or  motor  operating  conditions.  The  work  of  Kim 
[Ref.  9]  concentrated  on  the  development  of  the  diagnostic 
technique,  resulting  in  only  very  little  data  being  obtained. 

The  purposes  of  this  investigation  were  (1)  to  improve  the 
diagnostic  method  utilized  by  Kim  [Ref.  9]  by  significantly 
increasing  the  rate  of  acquisition  of  data  emd  by  increasing  the 
number  of  wavelengths  employed  and  (2)  to  apply  the  improved 
system  to  a  wider  range  of  test  conditions  to  determine  the 
effects  on  particle  size  and  optical  properties. 


II.  DZAGHOSTXC  METHOD 


A.  BACKCOtOUMD 

The  optical  analysis  of  particles  dates  back  to  Lord  Rayleigh 
and  the  nineteenth  century.  Rayleigh's  theory  explained  the 
scattering  of  light  from  particles  whose  dimensions  were  a  great 
deal  smaller  than  the  actual  wavelength  of  the  light.  Additionally 
this  allowed  him  to  explain  the  blue  sky  in  terms  of  the  small 
particles  suspended  in  the  upper  atmosphere.  (The  more  general 
technique  used  to  measure  all  particles  sizes  is  the  Mie  theory.) 
[Ref.  4,8] 

Since  Rayleigh's  initial  discoveries,  lasers  and  computers 
have  significantly  refined  the  optical  process,  but  the 
fundamental  approach  is  still  the  same.  A  collimated  light  beam 
is  sent  through  a  particle  cloud  and  information  about  the 
particles  is  determined  based  upon  the  wavelength  dependence  of 
the  scattering  and  extinction  of  the  light  beam  [Ref.  4]. 
Extinction  in  this  case,  refers  to  the  reduction  in  intensity  of 
the  light  due  to  absorption  and  scattering  of  the  incident  light 
by  the  particles. 

B.  THEORY 

The  Bouguer  or  Beer-Lambert  law  [Ref.  3,7],  describes  the 
exponential  decay  of  the  intensity  of  a  light  beeim  as  it  passes 
through  a  particle  field  of  uniform  size: 
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where  T  =  fraction  of  light  transmitted, 

Q  =  dimensionless  extinction  coefficient, 

A  =  cross  sectional  area  of  a  particle, 
n  =  number  concentration  of  particles, 

L  =  path  length, 

=mass  concentration  of  particles 
p  =  density  of  and  individual  particle 
d  =  particle  diameter 

Mie  scattering  theory  explains  the  interaction  of  a  plane 
electromagnetic  wave  with  a  spherical  object,  and  allows 
computation  of  Q  as  a  function  of  particle  size,  wavelength  of 
light,  and  complex  index  of  refraction  of  the  particle.  Since 
Mie's  initial  discoveries,  further  extensions  of  his  theory  have 
been  made,  making  it  applicable  to  nonspherical  particles  as  well. 
The  theory  holds  for  nonspherical  particles  if  they  are 
approximated  by  spherical  particles  whose  diameters  give  them  the 
same  volume  as  the  nonspherical  particles.  [Ref.  3,7] 

Dobbins  [discussed  in  Ref.  7]  revised  Bouguer's  transmission 
law  for  a  polydisperse  system  of  particles; 
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where  Q  is  the  average  extinction  coefficient,  and  d32  is  the 
voiume-to-surface  (Sauter)  mean  particle  diameter.  Taking  the 
natural  logarithm  of  equation  (2)  yields; 


2  P42 


(3) 


Since  the  light  beam  has  identical  path  lengths,  p,  C^,  and  d32  for 
each  wavelength  through  the  exhaust  plume,  the  ratio  of  the  In- 
transmissions  at  any  two  wavelengths  is  equal  to  the  ratio  of  the 
calculated  extinction  coefficients  for  the  saune  wavelengths; 


iAr<ly} 

A  computer  code  [Ref.  7],  utilizing  Mie  theory,  was  used  to 
generate  the  values  for  “q. 
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The  transmittances  at  each  of  the  various  wavelengths  are 
determined  experimentally.  The  procedure  requires  calculations  of 
o’  for  various  values  of  jL,  6^21  m(=n-ik),  and  an  assumed  particle 
size  distribution.  If  the  chosen  values  are  correct,  calculated 
extinction  coefficient  ratios  will  match  the  measured  In 
transmission  ratios  for  each  of  the  wavelengths  utilized.  Studies 
have  shown  little  variation  in  the  measured  indices  of 
refraction  with  wavelengths  in  the  visible  and  near  infrared 
ranges;  thus,  for  the  Mie  calculations,  the  indices  were  assumed 
constant  [Ref.  7]. 

Depending  upon  the  assumed  shape  of  the  size  distribution 
(monomodal,  bi-modal,  log  normal,  upper-limit, etc. ) ,  one  or  more 
parameters  are  required  to  characterize  it.  Based  upon  the 
studies  of  Chippett  and  Gray  [discussed  in  Ref.  7],  small 
particles  in  plumes  often  have  a  monomodal  log-normal 
distribution.  In  this  case  the  distribution  is  specified  by  d32 
and  the  standard  deviation (a) .  If  the  distribution  is  log-normal 
and  the  particles  are  absorbing  (complex  refractive  index)  then 
there  are  four  unknowns;  n,k,d32,  and  a,  requiring  a  minimum  of 
four  independent  measurements  involving  these  parameters. 

Calculations  with  the  Mie  code  showed  an  insensitivity  to 
the  absorption  index  for  the  reported  values  for  aluminum 
oxide (<10“^) ;  so  for  the  calculations,  the  absorption  index  of  the 
AI2O3  particles  was  assumed  to  be  zero  [Ref.  9].  With  three 
unknown  variables (d32,a,  and  n  (or  m) ) ,  a  minimum  of  three 
independent  In-transmission  ratios  were  needed.  The  actual 
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technique  utilizes  many  more  In-treuismission  ratios.  The  values 
of  n,  <332,  and  a,  are  found  by  regression  analysis,  i.e.  finding 
the  best  agreement  between  the  measured  and  calculated  results. 
Increasing  the  nximber  of  wavelengths,  increases  the  confidence 
level  in  the  data  correlation. 
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III.  EXPERIMENTAL  APPARATUS 


A.  ROCKET  MOTOR 

The  rocket  motor  used  in  this  experiment  was  the  Scime  one 
used  by  Kim  [Ref.  9].  The  external  dimensions  of  the  motor  were 
approximately  10.5”  in  length  by  3.25”  in  diameter.  The  motor 
consisted  of  an  ignitor,  nozzle,  burst  disk,  pressure  transducer, 
main  body  and  end  covers.  These  components  are  shown  in  figure  1. 


Figur*  1 

Motor  components 


The  Air  Force  Phillips  Laboratory  provided  a  GAP/AP  propellant 
containing  4.7%  alximinuro  for  use  in  the  motor.  Morton  Thiokol 
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provided  an  HTPB/AP  propellant  with  16%  altuninum.  The  propellants 
were  cut  to  burn  either  radially  (with  a  2”  outside  diameter,  a 
.375''web  and  a  length  of  1-2”)  or  axially (2”  outside  diameter  and 
a  length  of  1-1.5”).  The  portion  of  the  motor  walls  to  be  in 
contact  with  the  propellant  surface  were  lightly  coated  with  732- 
RTV  self-vulcanizing  silicone,  in  order  to  hold  the  propellant  in 
place,  and  to  act  as  an  burning  inhibitor.  To  ignite  the  motor, 
a  voltage  provided  by  a  12V  battery  was  applied  to  a  BKNO3 
ignitor . 

To  prevent  possible  overpressurization  of  the  motor  during 
firing,  a  disposable  1500  psi  biurst  disk  was  installed  on  the  side 
of  the  motor  wall.  To  record  motor  pressure,  a  pressure  port  for 
a  pressure  transducer  was  placed  in  the  combustion  chamber [Ref. 
9].  The  nozzle  specifications  are  shown  in  TABLE  1. 

B.  LIGHT  TRAM8MITTIMG  AMD  COLLECTZON  APPARATUS 

The  light  sovirce  for  the  experiment  consisted  of  an  Oriel 
Mercury-Xeon  white  light (Model  66085,  with  peak  wavelengths  at 
253.4,  302.2,  312.6,  334.2,  366.3,  404.7,  435.8,  546.1,  and 
577. Onm)  A  collimated  white  light  was  transmitted  through  the 
plume  and  onto  an  Oriel  Model  77400  Multispec  spectrograph  through 
a  25nm  wide  slit.  The  small  slit  width  was  necessary  to  insure 
that  the  amount  of  scattered  light  received  by  the  detector  was 
insignificant  in  comparison  to  the  total  amount  of  light  that 
penetrated  through  the  slit.  [Ref.  10]  The  light  was  then 
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TABLE  1  NOZZLE  SPECIFICATIONS 


reflected  from  a  mirror  to  a  diffraction  grating.  The  light  was 
then  reflected  off  of  another  mirror  to  a  1024  element  diode  array 
[Ref.  11].  A  schematic  of  the  optical  pattern  is  shown  in  Figure 

2. 


Figure  2 

Spectrograph  Optical  Pattern 
[Ref,  10] 


The  response  characteristic  of  the  photodiode  area  is  shown  in  fig 
3.  With  these  response  characteristics,  six  of  the  ten  available 
wavelengths  could  be  detected  (334.2,  366.3,  404.7,  435.8,  546.1, 
577. Onm) [Ref .11] . 


12 


RELATIVE  RESPONSE 


Figur*  3 

Diode  Response  Char. 
[Ref.  11] 


IV.  DATA  ACQUISITICaf  AMD  REDDCTION 

The  output  of  the  diode  array  was  sent  to  a  HP6942A 
Multiprogrammer,  where  the  analog  signal  was  converted  to  digital, 
and  stored  on  high  speed  memory  ceurds.  Collection  and  storage  of 
the  data  was  controlled  by  using  a  HP9836S  Computer.  The  timing 
circuit  of  the  diode  array  board  was  synchronized  with  the 
Multiprogrammer.  The  Multiprogrammer  system  consisted  of: 

1  SOOKhz  Analog  to  Digital  Converter  Card 

1  High  Speed  Memory  Card 

4  Memory  Expansion  Cards 

1  Timer  Pacer  Card 

The  A/D  card  was  configured  for  external  triggering, 
allowing  the  timer  pacer  card  to  control  the  rate  of  conversion. 
The  timer  pacer  card  was  also  configiured  for  external  triggering, 
allowing  it  to  be  paced  by  the  blanking  pulse  from  the  diode  array 
board.  With  this  control,  data  were  collected  at  the  beginning  of 
each  sweep  of  the  array,  without  large  gaps  in  data  between 
sweeps.  A  timing  diagram  is  shown  in  figure  4. 

The  data  from  the  sweeps  were  stored  in  bineury  data 
files(BDAT).  Each  BDAT  contained  data  from  eight  sweeps  of  the 
1024  element  diode  array(809€  bytes) .  There  were  six  peak 
voltages  in  each  sweep,  corresponding  to  each  of  the  wavelengths. 
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During  a  single  test  four  sequential  BOAT  files  were  used,  giving 
a  total  of  32  sweeps  of  the  diode  array.  Several  progreuns  were 
written  [and  adapted  from  Ref.  12] to  aid  in  finding  the  peaks  and 
to  determine  if  emy  of  the  sweeps  contained  bad  data.  The  peak 
voltages  were  averaged  according  to  wavelength  for  both  p2urticle 
and  no  particle  data.  The  tremsmittances  were  determined  by 
taking  the  ratio  of  particle  to  no  peurticle  data  for  each 
wavelength. 


Clock 
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Figure  4 

Timing  Diagram 
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y.  RESULTS  AND  DISCUSSION 


A  summary  of  the  propellants,  nozzle  sizes,  operating 
pressures,  and  data  collection  locations  for  each  run  is  shown  in 
TABLE  2.  Various  nozzle  geometries  were  used  to  provide  different 
chcunber  pressures  and  different  degrees  of  over-and  under¬ 
expansion. 

A.  DATA  REDUCTION 

The  index  of  refraction  of  the  AI2O3  peurticles  in  the  plume  was 
assumed  to  be  independent  of  wavelength,  and  the  absorption  index 
was  taken  to  be  zero  for  the  reasons  discussed  above.  Q  was 
determined  as  a  function  of  d32,  n,  and  a  using  the  Mie  code.  The 
ratio  of  mean  extinction  coefficients  were  then  plotted  against 
the  In-transmission  ratios  for  each  pair  of  wavelengths,  for  a 
total  of  fifteen  ratios.  If  the  ratios  had  corresponded  with  each 
other  perfectly,  the  data  points  would  have  formed  a  forty-five 
degree  line  passing  through  the  origin.  Short  of  this,  the  best 
fit  of  data  would  have  the  least  scatter  about  this  line,  and 
would  form  a  line  with  a  slope  approximately  equal  to  one. 

Finding  the  correct  combination  of  d32,  a,  and  n  was  an 
iterative  process  that  beg£Ui  with  expected  values  taken  from  Kim's 
report  [Ref.  9].  A  linear  least-squares  fit  method  was  used  to 
find  the  correlation  coefficient,  R,  as  a  measure  of  the 
relationship  between  the  In-transmission  ratios  and  the  extinction 
coefficient  ratios.  The  square  of  the  correlation  coefficient 
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expressed  as  a  percentage  was  used  to  measure  how  closely  the 
data  points  matched  the  regression  line.  The  value  closest  to 
one  indicated  the  best  linear  correlation  between  the  data.  The 
X-Coeff icient  (slope)  of  the  regression  analysis,  indicated  how 
closely  the  slope  of  the  least-squares  fit  matched  a  45° 
line. [Ref.  13:p.  17] 


-17 


TABLE  2.  TESTS  CX)NDUCTED 


Test 

# 

%  AL 

GC 

GL 

(in.) 

dth 
(in. ) 

Pc  . 
(psia) 

Pe 

(psia) 

PL 

x(x/de) 

r(r/re) 

(in.) 

1 

16% 

rad 

2 

.265 

.811 

780 

10.9 

slightly 

over- 

expanded 

6(7.4) 
1.4(3. 5) 

2 

16% 

rad 

1.5 

.265 

.811 

405 

5.5 

over- 

expanded 

6. 5(8.0) 
1.4(3. 5) 

3 

4.69% 

end 

1 

.217 

.381 

630 

42.3 

under¬ 

expanded 

6(15.7) 
1.2(6. 3) 

4 

4.69% 

end 

1 

.265 

.811 

195 

2.7 

highly 

over- 

expanded 

6(7.4) 

1. 2(3.0) 

end  -  end  burning  grain 

GC  -  grain  configuration 

GL  -  grain  length 

rad  -  radial  burning  grain 

PL  -  plume  location  aft  of  nozzle  exit 
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B.  MEASDREMEMT  ACCURACY 


Earlier  studies  conducted  by  Kim  [Ref.  9]  using  the  multiple 
wavelength  transmittance  apparatus  had  determined  the  approximate 
experimental  measurement  accuracy  of  the  transmittances  (due  to 
electronic  noise)  to  be  within  0.5%.  When  wavelengths  are  closely 
spaced  and  or  higher  transmittances  (above  90%)  occur,  this 
accuracy  results  in  a  large  uncertainty  in  the  In-transmission 
ratio.  To  reduce  the  effects  of  this  problem  and  increase  the 
confidence  in  the  measurement  accuracy,  the  number  of  wavelength 
ratios  was  increased.  The  larger  number  of  closely  spaced  data 
points  reduce  the  possibility  of  fitting  several  regression  lines 
through  the  same  data.  In  addition,  the  higher  speed  data 
acquisition  permitted  32  sweeps  of  the  diode  array  in  70msec. 
This  significantly  increased  the  number  of  repeated  measurements 
for  averaging  and  decreased  the  total  sampling  time  from  that 
utilized  by  Kim  [Ref.  9].  The  latter  results  in  less  variation 
due  to  changing  motor  conditions  during  the  measurement  sequence. 
When  the  best  solution  was  found  using  fifteen  wavelength  ratios, 
d32,  a,  and  n  were  perturbed  about  their  nominal  values  until  a  2% 
variation  in  slope  (X-coefficient)  was  noted  in  the  plotted 
solution.  The  graphical  results  of  these  perturbations  are  shown 
in  Figvures  6  through  17.  From  this  analysis,  for  a  2%  variation 
in  X-Coefficient,  the  accuracy  of  the  measured  values  of  d32,  n, 
and  a  were  estimated.  A  summary  of  these  results  is  shown  in 
Table  3.  The  resulting  average  uncertainty  in  d32  and  a  were 
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slightly  increased  from  that  determined  by  Kim.  However,  the 
uncertainty  in  the  index  of  refraction  was  significantly  reduced. 
The  good  correlation  of  the  data  indicated  that  the  particle  size 
distribution  could  be  satisfactorily  represented  by  a  log-normal 
distribution. 
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Figure  6. 

Results  of  Perturbation  of  <132  about  .33, Run  1 
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Figure  7 . 

Results  of  Perturbation  of  <132  about  .21, Run  2 
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Figure  9. 

Results  of  Perturbation  of  d32  about  .33, Run  4 
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Figure  10. 

Results  of  Perturbation  of  n  about  1.63, Run  1 
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Figure  11. 

Results  of  Perturbation  of  n  about  1.63, Run  2 
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Figure  12. 

Results  of  Perturbation  of  n  about  1.67, Run  3 
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Figure  13. 

Results  of  Perturbation  of  n  about  1.61, Run  4 
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Figure  14. 

Results  of  Perturbation  of  sigma  about  1.4, Run  1 
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Figure  15. 

Results  of  Perturbation  of  sigma  about  1.47, Run  2 
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Figure  16. 

Results  of  Perturbation  of  sigma  about  1.6, Run  3 
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Figure  17 . 

Results  of  Perturbation  of  sigma  about  1.6, Run  4 
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TABLE  3.  TEST  RESULTS 


Test 

a 

n 

r2 

X-Coef 

1 

.33+.  02/1 

1.40+.10 

1.63+.02 

.99381 

.999540 

2 

.21+.  02/1 

.98659 

1.00479 

3 

.33+. 02m 

1.60+.20 

1.67+.03 

.99676 

.995403 

4 

BBHB!!! 

.97354 

1.00699 

aver 

uncertainty 

.30+.  02/1 

1.52+.12 

1.64+.04 

C.  EFFECTS  OF  TEST  CONDITIONS 

For  the  16%  AL  propellant,  the  change  in  chamber  pressure  and 
degree  of  over-expansion  did  not  have  a  profound  effect  on  the 
optical  properties.  The  index  of  refraction  remained  the  same, 
and  the  change  in  the  standard  deviation  of  the  log-normal 
distribution  was  within  the  measurement  uncertainty.  The  Sauter 
mean  diameter  decreased  in  run  2;  however,  the  measurement  was 
made  further  aft  in  the  plume.  D32  can  change  rapidly  with  plume 
position  [Ref.  14]. 

Significantly  different  pressures  and  expansion  ratios  had 
little  effect  on  the  properties  \dien  using  the  4.69%  AL.  The  16% 
AL  propellant  contained  a  small  amount  of  iron  oxide  as  a  burn 
rate  catalyst  and  used  an  HTPB  binder.  The  4.69%  AL  propellant 
had  no  iron  oxide  and  used  a  GAP  binder.  With  these  significant 
differences,  the  value  of  the  index  of  refraction  remained  at 
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approximately  1.64,  essentially  the  same  as  measured  by  Kim  [Re? 
9].  This  result  indicates  that  the  small  particles  that  are 
rapidly  cooled  in  the  nozzle  flow  are  probably  Y-AI2O3  as  observed 


by  other  investigators  [Ref.  15,16,17].  a-Al203  has  a 

significantly  higher  index  of  refraction. 


34 


VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  increase  in  the  number  of  wavelengths  (5  to  6)  utilized  by 
the  multiple  wavelength  transmittance  apparatus  provided  more  In- 
transmittance  ratios  and  increased  confidence  in  the  measured 
parameters  obtained  by  least-squares  fitting  of  the  data.  The 
increased  A/D  rate  (33KHz  to  SOOKHz)  resulted  in  the  ability  to 
obtain  more  sweeps  (32  vs  8)  in  significantly  less  time(70ms  vs 
250ms) .  This  resulted  in  two  benefits.  First,  time  variations  in 
motor  operating  conditions  had  much  less  effect.  Second,  more 
data  were  obtained  for  averaging,  thereby  reducing  the  effects  of 
the  elimination  of  obviously  bad  data  (due  to  large  particles 
being  present  during  one  sweep,  etc.) .  The  rew  computer  code  also 
eliminated  the  need  for  hand  selection  of  the  peeik  voltages  from 
the  diode  array,  significantly  reducing  data  reduction  time. 
These  modifications  did  not  improve  the  reported  accuracy  for 
determining  <1^2  and  a  [Ref.  9],  but  did  significantly  improve  the 
accuracy  for  determining  the  index  of  refraction  (8%  to  3%) .  The 
average  values  of  ^22,  n,  and  a  obtained  for  the  two  aluminized 
propellants  were  d32=.30+.02^,  n  =  1.64+.04  and  a  =  1.52+.12. 

The  results  of  this  investigation,  together  with  that  of 
Ref. 9,  indicate  that  the  small  aluminum  oxide  particles  (that 
dominate  the  outer  plmne  regions)  are  ^“^1203#  independent  of 
propellant  composition,  motor  operating  conditions  and  nozzle 
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geometry.  In  addition,  the  good  correlation  of  the  data  indicated 
that  the  small  particles  can  be  adequately  represented  by  a 
monomodal,  log-normal  distribution. 

Subsequent  work  should  consider  further  modification  of  the 
Fortran  code  used  the  generate  the  Mie  extinction  coefficients. 
Currently  the  code  generates  the  extinction  coefficient  ratios  2md 
d32  for  three  wavelengths.  Increasing  the  number  of  wavelengths 
to  six,  and  adding  a  regression  function  to  automatically 
determine  the  best  correlation  between  the  measured  and  calculated 
results,  would  shorten  data  reduction  time  and  increase  the 
overall  accuracy. 
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APPENDIX  A  -  SYSTEM  UPGRADES 


The  upgrades  to  the  system  used  by  Kim  [Ref.  9], 
included  increasing  the  speed  of  A/D  conversion  from  33kh2  to 
SOOkhz.  This  involved  replacing  the  HP69751A  A/D  card,  with  the 
HP69759A  A/D  card,  and  adjusting  the  HP69736A  Timer/Pacer  card 
accordingly.  Additionally,  the  HP69790B  memory  cards  were 
replaced  with  a  HP69791A  High  Speed  Memory  card,  and  four 
HP69792A  Memory  Expansion  cards.  This  exchange  increased  the 
system  storage  capacity  from  less  than  14  thousand  words  to  over 
850  thousand  words  of  high  speed  memory.  In  essence,  the  new 
system  permitted  the  collection  of  sixty  times  as  much  data  in  one 
fifth  the  time. 

The  upgrade  in  computer  hardware  also  necessitated 
updating  the  software  programs  used  to  control  data  collection, 
storage,  and  reduction.  The  original  system  stored  data  to  one  of 
two  5  1/4  floppy  drives,  up  to  the  capacity  of  the  diskettes.  The 
increased  storage  capacity  of  the  memory  cards  necessitated  moving 
permanent  data  storage  to  the  hard  drive  of  the  HP9836S  computer. 
Currently,  the  data  collection  and  storage  program  saves  eight 
binary  data  files  of  data  at  a  time.  The  procedure  used  for 
creating  and  storing  the  files  is  readily  expandable  to  15-20 
files  of  12  sweeps  each.  The  slow  transfer  rate  of  the  HP9836S 
computer,  both  from  the  memory  cards  to  the  computer,  and  from  the 
computer  to  the  hard  drive,  is  why  only  eight  files  are  currently 
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being  created  and  transferred.  An  HP  interface  card  is  being 
purchased  for  a  new  486PC.  This  will  permit  the  full  potential  of 
the  newer  memory  system  to  be  realized.  The  revised  data 
collection  program  is  shown  in  Appendix  C. 
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APPENDIX  B  -  HEWLETT  PACKARD  lfDLTIPR06RAMMER(HP6942A) 


1.  GENERAL  DESCRIPTION 

The  HP  6942A  Multlprogrammer  system  is  a  self-contained 
module  that  can  be  used  in  a  wide  variety  of  data  capture 
applications [Ref.  18:p.  l-l,Ref.  19:p.  32].  The  Multiprogrammer 
utilizes  a  number  of  I/O  cards,  and  is  designed  to  operate  on  the 
HP  Interface  Bus  under  control  of  a  desk  top  computer.  The  main 
advantage  of  the  Multiprogrammer  and  its  I/O  Ceurds  over  other  data 
acquisition  systems  is  that  its  cards  can  be  controlled  by  dual 
means  [Ref.  18 :p.  3-1].  A  I/O  card  can  be  addressed  cuid 
controlled  by  the  computer  via  the  back  plane  edge  of  the  card,  or 
external  events  can  control  operation  of  the  ceurd  via  the  edge 
connectors  on  the  cards.  Both  means  of  controlling  the  I/O  ceurds 
to  collect  data  have  been  used.  The  A/D  card  and  the  Timer/Pacer 
card  were  both  triggered  externally  by  inputs  from  the  photodiode; 
however  the  actual  data  acquisition  process  was  initiated  by  the 
computer  via  inputs  from  the  motor  pressiure  treuisducer  voltage. 

The  Multiprogrammer  can  be  configured  to  operate  with  of  up  to 
sixteen  I/O  cards,  depending  upon  its  intended  purpose.  The 
current  arrangement  emphasized  the  capturing  of  analog  data  for 
conversion  and  storage  to  the  heurd  drive  of  the  HP9836S 
computer. [Ref .  18 :p.  2-1] 
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2.  500  KHZ  A/D  CARD  (HP69759A) 

As  shipped  from  the  factory,  the  Analog  to  Digital  conversion 
card  had  a  nominal  +  10-volt  unity-gain  input  range  with  5- 
millivolt  resolution.  It  can  be  cycled  either  internally  (via  the 
back  plane)  or  externally  (via  the  edge  connector) .  When  it  is 
triggered,  the  A/D  card  converts  the  input  voltage  into  a  12-bit, 
two's  complement  binary  value  within  two  microseconds.  [Ref.  20;p. 
1-2]  The  card  has  14  programmedsle  mode  functions  which  include, 
gain  , trigger,  and  external  lockout  control.  In  this 
investigation,  the  card  was  operated  in  the  +  1  range,  received 
input  voltages  from  the  diode  array,  auid  was  triggered  externally 
by  the  timer /pacer  card. 

3.  RIOR  SPEED  MEMORY  (69791A)  AMD  MEMQRY  EZPAMSION  CARDS  (HP69792A) 

The  Model  69791A  High  Speed  Memory  card  consists  of  two 
printed  circuit  boards  that  are  physically  attached  to  each  other 
by  four  teflon  snap  connectors.  The  card  has  dynamic  Random 
Access  Memory  (RAM)  capable  of  storing  up  to  65,536  16-bit  words  at 
high  speed  [Ref.  21 :p.  1-1].  The  16-bit  data  can  be  written  to, 
or  read  from,  the  card  by  either  the  computer  or  an  external 
device.  Because  the  card  is  capable  of  simultaneously  reading  and 
writing  data,  the  computer  can  transfer  data  via  the  back  plane  of 
the  Multiprogrammer,  vdiile  external  devices  transfer  data  via  the 
edge  connector. 

The  Model  69792A  Memory  Expansion  card  is  electrically 
connected  to  the  HP  69791A  via  a  chaining  cable.  Each  Memory 
Expansion  card  that  is  added  to  the  Memory  Ceurd  System  increases 
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the  number  of  16-bit  words  that  can  be  stored  by  196,608  [Ref. 
21 :p.  1-1].  Up  to  five  expansion  cards  can  be  used  with  one  HP 
69791A  for  a  storage  capacity  of  1,048,576  words.  The  present 
configuration  uses  four  expansion  cards,  for  a  memory  capacity  of 
851,968  words.  The  Memory  Card  System  has  two  basic  modes  of 
operation:  FIFO  (First  In/First  Out)  mode  and  Recirculate  mode. 
In  the  current  setup  the  Memory  Card  System  was  configured  in  the 
FIFO  mode,  and  acted  as  a  data  buffer.  Converted  voltages  from 
the  A/D  card  were  stored  in  the  Memory  Card  System  before  and 
during  the  actual  rocket  motor  firing,  and  later  transferred  to 
the  hard  drive  of  the  HP9836S  cou^uter. 

4.  TINER/PACER  CARD(HP69736A) 

The  Model  HP  697 3 6A  Timer /Pacer  card  is  a  square  wave  pulse 
generator  that  can  generate  output  pulses  from  1  ^sec  to  18.2  hrs 
in  duration  [Ref.  22 :p.  1-1],  The  card  operates  in  either  the 
one-shot  mode  (single  output  pulse)  or  the  recirculate  mode 
(continuous  train  of  output  pulses) .  In  the  present  conf  igvuration 
the  Timer/Pacer  was  operated  in  the  recirculate  mode  and  was 
triggered  externally  by  the  blemking  pulse  of  the  diode  array. 
This  synchronized  the  card  pulse  train  with  the  first  sweeps  of 
the  diode  array,  which  in  turn  controlled  the  conversion  rate  of 
data  by  the  A/D  card. 
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